Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe.

BACKGROUND: Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry. METHODS AND RESULTS: MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank. CONCLUSIONS: Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Positive selection on two mitochondrial coding genes and adaptation signals in hares (genus Lepus) from China.

BACKGROUND: Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection. RESULTS: Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species. CONCLUSIONS: Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient.

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Genetic diversity of the toll-like receptor 2 (TLR2) in hare (Lepus capensis) populations from Tunisia.

Toll-like receptors (TLRs) are a major group of proteins that recognize molecular components of infectious agents, known as pathogen associated molecular patterns (PAMPs). The structure of these genes is similar and characterized by the presence of an ectodomain, a signal transmembrane segment and a highly conserved cytoplasmic domain. The latter domain is homologous to the human interleukin-1 receptor (IL1R) and human IL-18 receptor (IL-18R) and designated TIR domain. The latter domain of the TLR genes was suggested to be very conservative and its evolution is driven by purifying selection. Variability and evolution of the TIR sequences of TLR2 gene were studied in three hare populations from Tunisia with different ecological characteristics (NT-North Tunisia with Mediterranean, CT-Central Tunisia with semi-arid, and ST-South Tunisia with arid climate). Sequencing of a 372bp fragment of TIR2 revealed 25 alleles among 110 hares. Twenty variable nucleotide positions were detected, of which 7 were non-synonymous. The highest variability was observed in CT, with 16 polymorphic positions. In ST, only 4 polymorphic nucleotide positions were detected with all diversity values lower than those recorded for the other two populations. By using several approaches, no positive selection was detected. However, evidence of purifying selection was found at two positions. The logistic models of the most common TIR2 protein variant that we run to examine whether its occurrence was affected by climatic variation independent of the geographic sample location suggested only a longitudinal effect. Finally, the mapping of the non-synonymous mutations to the inferred tertiary protein structure showed that they were all localized in the different loop regions. Among all non-synonymous substitutions, three were suggested to be deleterious as evidenced by PROVEAN analysis. The observed patterns of variability characterized by low genetic diversity in ST might suggest that the TIR region was more affected, than other markers, by genetic drift or/and that these patterns were shaped by different selective pressures under different ecological conditions. Notably, this low diversity was not detected by other (putatively neutral) microsatellite markers analysed in the course of other studies. But low diversity was also found for two MHC class II adaptive immune genes. As expected from functionally important regions, the evolution of the TIR2 domain is mainly driven by purifying selection. However, the occurrence of deleterious non-synonymous substitutions might highlight the flexible evolution of the TIR genes and/or their interactions with other proteins.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia.

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

Selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes in hares (Lepus capensis L., 1758) from a steep ecological gradient in North Africa.

BACKGROUND: Recent studies of selection on mitochondrial (mt) OXPHOS genes suggest adaptation due mainly to environmental variation. In this context, Tunisian hares that display several external phenotypes with phylogenetically rather homogenous gene pool and shallow population structure provide a good precondition to detect positive selection on mt genes related to environmental/climatic variation, specifically ambient temperature and precipitation. RESULTS: We used codon-based methods along with population genetic data to test for positive selection on ATPase synthase 6 (ATP6) and NADH dehydrogenase 2 (ND2) of cape hares (Lepus capensis) collected along a steep ecological gradient in Tunisia. We found significantly higher differentiation at the ATP6 locus across Tunisia, with sub-humid Mediterranean, semi-arid, and arid Sahara climate than for fourteen unlinked supposedly neutrally evolving nuclear microsatellites and mt control region sequences. This suggested positive selection on ATP6 sequences, which was confirmed by several codon-based tests for one sequence site that together with a second site translated into four different amino acids. Positive selection on ND2 sequences was also confirmed by several codon-based tests. The corresponding frequencies of the two most prevalent variants at each locus varied significantly across climate regions, and our logistic general linear models of occurrence of those proteins indicated significant effects of mean annual temperature for ATP6 and mean minimum temperature of the coldest month of the year for ND2, independent of geographical location, annual precipitation, and the respective co-occurring protein at the second locus. Moreover, presence of the ancestral ATP6 protein, as inferred from phylogenetic networks, was positively affected by the simultaneous presence of the derived ND2 protein and vice versa, independent of temperature, precipitation, or geographic location. Finally, we obtained a significant coevolution signal for the ancestral ATP6 and derived ND2 sequences and vice versa. CONCLUSIONS: positive selection was strongly suggested by the population genetic approach and the codon-based tests in both mtDNA genes. Moreover, the two most prevalent proteins at the ATP6 locus were distributed at significantly varying frequencies across the study area with a significant effect of mean annual temperature on the occurrence of the ATP6 proteins independent of geographical coordinates and the co-occuring ND2 protein variant. For ND2, occurrence of the two most frequent protein variants was significantly influenced by the mean minimum temperature of the coldest month, independent of the co-occurring ATP6 protein variant and geographical coordinates. This strongly suggests direct involvement of ambient temperature in the adaptation of the studied mtOXPHOS genes.

Variation of partial transferrin sequences and phylogenetic relationships among hares (Lepus capensis, Lagomorpha) from Tunisia.

North African hares are currently included in cape hares, Lepus capensis sensu lato, a taxon that may be considered a superspecies or a complex of closely related species. The existing molecular data, however, are not unequivocal, with mtDNA control region sequences suggesting a separate species status and nuclear loci (allozymes, microsatellites) revealing conspecificity of L. capensis and L. europaeus. Here, we study sequence variation in the intron 6 (468 bp) of the transferrin nuclear gene, of 105 hares with different coat colour from different regions in Tunisia with respect to genetic diversity and differentiation, as well as their phylogenetic status. Forty-six haplotypes (alleles) were revealed and compared phylogenetically to all available TF haplotypes of various Lepus species retrieved from GenBank. Maximum Likelihood, neighbor joining and median joining network analyses concordantly grouped all currently obtained haplotypes together with haplotypes belonging to six different Chinese hare species and the African scrub hare L. saxatilis. Moreover, two Tunisian haploypes were shared with L. capensis, L timidus, L. sinensis, L. yarkandensis, and L. hainanus from China. These results indicated the evolutionary complexity of the genus Lepus with the mixing of nuclear gene haplotypes resulting from introgressive hybridization or/and shared ancestral polymorphism. We report the presence of shared ancestral polymorphism between North African and Chinese hares. This has not been detected earlier in the mtDNA sequences of the same individuals. Genetic diversity of the TF sequences from the Tunisian populations was relatively high compared to other hare populations. However, genetic differentiation and gene flow analyses (AMOVA, FST, Nm) indicated little divergence with the absence of geographically meaningful phylogroups and lack of clustering with coat colour types. These results confirm the presence of a single hare species in Tunisia, but a sound inference on its phylogenetic position would require additional nuclear markers and numerous geographically meaningful samples from Africa and Eurasia.

Sex-specific selection for MHC variability in Alpine chamois.

BACKGROUND: In mammals, males typically have shorter lives than females. This difference is thought to be due to behavioural traits which enhance competitive abilities, and hence male reproductive success, but impair survival. Furthermore, in many species males usually show higher parasite burden than females. Consequently, the intensity of selection for genetic factors which reduce susceptibility to pathogens may differ between sexes. High variability at the major histocompatibility complex (MHC) genes is believed to be advantageous for detecting and combating the range of infectious agents present in the environment. Increased heterozygosity at these immune genes is expected to be important for individual longevity. However, whether males in natural populations benefit more from MHC heterozygosity than females has rarely been investigated. We investigated this question in a long-term study of free-living Alpine chamois (Rupicapra rupicapra), a polygynous mountain ungulate. RESULTS: Here we show that male chamois survive significantly (P = 0.022) longer if heterozygous at the MHC class II DRB locus, whereas females do not. Improved survival of males was not a result of heterozygote advantage per se, as background heterozygosity (estimated across twelve microsatellite loci) did not change significantly with age. Furthermore, reproductively active males depleted their body fat reserves earlier than females leading to significantly impaired survival rates in this sex (P < 0.008). This sex-difference was even more pronounced in areas affected by scabies, a severe parasitosis, as reproductively active males were less likely to survive than females. However, we did not find evidence for a survival advantage associated with specific MHC alleles in areas affected by scabies. CONCLUSIONS: Increased MHC class II DRB heterozygosity with age in males, suggests that MHC heterozygous males survive longer than homozygotes. Reproductively active males appear to be less likely to survive than females most likely because of the energetic challenge of the winter rut, accompanied by earlier depletion of their body fat stores, and a generally higher parasite burden. This scenario renders the MHC-mediated immune response more important for males than for females, which implies a relatively stronger selection pressure on MHC genes in males than in females.

Mitochondrial CR-1 variation in Sardinian hares and its relationships with other Old World hares (Genus Lepus).

Among the European fauna, the Sardinian hare (Lepus sp.) is peculiar in that it differs from all other hares inhabiting the continent. Here, we report on the variation of a 461 bp sequence of hypervariable domain 1 of the mitochondrial control region, examined in 42 hares collected throughout Sardinia and compared to the corresponding sequences of different Lepus taxa. Seventeen novel haplotypes were found in the Sardinian population, resulting in a haplotype diversity of 0.840 and a nucleotide diversity of 0.012. As a result of Bayesian and principal coordinates analyses, Sardinian hares were grouped with North African hares, constituting a monophyletic clade that diverges from all other Old World hares, including Cape hares from South Africa and East Asia. Hence, our data agree that populations inhabiting North Africa and Sardinia form a distinct taxon, which could possibly be included in the L. capensis superspecies. Moreover, two corresponding lineages can be found in Sardinia and Tunisia, providing evidence of a common origin of the two populations and thus supporting the hypothesis that North African hares were introduced into the island in historical times. Our data show that the two lineages differ in their geographic distribution throughout the island and that the wild Sardinian population also shows the signature of a postintroduction demographic expansion.

Biochemical genetic relationships among Tunisian hares (Lepus sp.), South African cape hares (L. capensis), and European brown hares (L. europaeus).

Tunisian hares (n = 45), currently assigned to Lepus capensis, were assayed for allelic variation at 40 allozyme loci, and allele frequencies at 32 loci were directly compared with earlier data of South African cape hares (L. capensis, n = 9) and European brown hares (L. europaeus, n = 244) to reveal genetic relationships among them. European mountain hares (L. timidus, n = 200) were used for outgroup comparison. In the Tunisian hares 27.5% of the loci were polymorphic with 2-4 alleles. Among all alleles at polymorphic loci, 15.1% occurred exclusively in Tunisian hares, 5.7% exclusively in cape hares, and 7.5% exclusively in brown hares at low frequencies. Not a single locus showed alternately fixed alleles between the samples of the L. capensis/L. europaeus complex. Levels of absolute and relative genetic differentiation among the samples of the L. capensis/L. europaeus complex were low, relative to pairwise comparisons involving mountain hares. Diverse cluster analyses and multidimensional scaling of various pairwise genetic distance matrices concordantly grouped Tunisian hares with brown hares, and South African cape hares clustered only slightly farther apart, whereas mountain hares were distinctly separate. These results suggest regionally distinct phylogenetic units within an overall cohesive gene pool in the L. capensis/L. europaeus complex, supporting Petter's view that all North African hares belong to L. capensis except for one local population of savanna hares, and that cape hares and brown hares are conspecific.